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q . Abstract 

•l-H . 

The characterization of non-stationary signals requires joint time 
and frequency information. However, time (t) and frequency (w) being 
Qh| non-commuting variables, there cannot be a joint probability density 

in the (t, to) plane and the time- frequency distributions, that have been 
proposed, have difficult interpretation problems arising from negative 
or complex values and spurious components. As an alternative we pro- 
ps) ' pose to obtain time-frequency information by looking at the marginal 

distributions along rotated directions in the (t,oj) plane. The rigor- 
ous probability interpretation of the marginal distributions avoids all 
t^- ' interpretation ambiguities. Applications to signal analysis and signal 

0\ ■ detection are discussed as well as an extension of the method to other 



> 
(N 



^ , pairs of non-commuting variables. 

•i-H 

43 ! 1 Introduction 

Non-stationary signals have a time-dependent spectral content, therefore, an 
adequate characterization of these signals requires joint time and frequency 
information. Among the many time-frequency (quasi) distributions [|T| |2| that 
have been proposed, Wigner-Ville's (WV)@ (§] 

W(t,u) = J f(t + |) /' (t - |) e— du (1) 
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for an analytic signal / (£), is considered to be optimal in the sense that it 
satisfies the marginals, it is time-frequency shift invariant and it possesses 
the least amount of spread in the time-frequency plane. 

However, the WV distribution has, in general, positive and negative val- 
ues and may be non-zero in regions of the time-frequency plane where either 
the signal or its Fourier transform vanish. Therefore, despite the fact that the 
WV distribution is an accurate mathematical characterization of the signal, 
in the sense that it can be inverted by 



its interpretation for signal detection and recognition is no easy matter, be- 
cause of the negative and " spurious" components. The origin of this problem 
lies in the fact that t and u being non-commuting variables, they cannot be 
simultaneously specified with absolute accuracy and, as a result, there cannot 
be a joint probability density in the time-frequency plane. Therefore no joint 
distribution, even if positive |5|], may be interpreted as a probability density. 

Looking back at the original motivation leading to the construction of the 
time-frequency distributions, namely the characterization of non-stationary 
signals, we notice that we are asking for more than we really need. To charac- 
terize a non- stationary signal what we need is time and frequency-dependent 
information, not necessarily a joint probability density, a mathematical im- 
possibility for non-commuting variables. The solution is very simple. The 
time density \f(t) | 2 projects the signal intensity on the time axis and the spec- 
tral density \f(u)\ 2 projects on the frequency axis. To obtain the required 
time-frequency information, all we need is a family of time and frequency 
functions sg(t,u), depending on a parameter £, which interpolates between 
time and frequency. Projecting the signal intensity on this variable, that is, 
computing the density along the s^— axis, one obtains a function 



that has, for each £, a probability interpretation. The simplest choice for 
is a linear combination 



the parameter £ being the pair (/i, v). For definiteness we may choose 





M( a ,fl = l/(*)| 2 



(3) 



S = lit + VUJ 



(4) 



T, f2 being a reference time and a reference frequency adapted to the signal 
to be studied. The function M(s,8) interpolates between \ f{t)\ 2 and \ f(u)\ 2 
and, as we will prove below, contains a complete description of the signal. 
For each 9 the function M(s, 9) is strictly positive and being a bona-fide 
probability (in s) causes no interpretation ambiguities. A similar approach 
has already been suggested for quantum optics|| and quantum mechanics J?]] 
H] II' the non-commuting variable pairs being respectively the quadrature 
phases (a r ,Oj) and the position-momentum (q,p). 

This approach, in which to reconstruct an object, be it a signal in signal 
processing or a wave function in quantum mechanics, one looks at its prob- 
ability projections on a family of rotated axis, is similar to the computerized 
axial tomography (CAT) method. The basic difference is that in CAT scans 
one deals with a pair (x, y) of commuting position variables and here we deal 
with a plane defined by a pair of non-commuting variables. For this reason 
we call the present approach non- commutative tomography (NCT). 

The paper is organized as follows. In Section 2 we construct the NCT 
signal transform and show its positivity and normalization properties. We 
also establish the invertibility of the transformation, which shows that it 
contains a complete description of the signal and establish its relation to 
the WV distribution. Because the NCT transform involves the square of 
the absolute value of a linear functional of the signal, it is actually easier to 
compute than bilinear transforms like WV. 

In Section 3 we work out the analytical form of the NCT transform for 
some signals and also display the M(s, 9) in some examples. We also deal with 
the problem of using NCT to detect the presence of signals in noise for small 
signal to noise ratios (SNR). Here the essential observation is that, for small 
SNR, the signal may be difficult to detect along t or lj, however, it is probable 
that there are other directions on the (t, u) plane along which detection might 
be easier. It is the consistent occurrence of many such directions that supplies 
the detection signature. 

Finally in Section 4 we point out that the NCT approach may also be used 
for other pairs of non-commuting variables of importance in signal processing. 
As an example we work out the relevant formulas for the scale-frequency pair. 
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2 Non- commutative time-frequency tomogra- 
phy 

Because the Fourier transform of a characteristic function is a probability 
density, we compute the marginal distribution for the variable s = fit + vlo 
using the characteristic function method. Frequency and time are operators 
acting in the Hilbert space of analytic signals and, in the time-representation, 
the frequency operator is uo = —id/dt . The characteristic function C(k) is 

C(ife) = ( e *0i*+^)) = J f*(t) e ^t-iud/dt) dt ^ 
where f(t) is a normalized signal 

/ \f(t)\ 2 dt = l 

The Fourier transform of the characteristic function is a probability density 

M (s, fi,u) = ^- I C{k)e~ iks dk (7) 

2 71 J 

After some algebra, one obtains the marginal distribution (0) in terms of the 
analytical signal 

2 

(8) 

with normalization 

J M(s,(i,u) ds = 1 (9) 

For the case \i = 1, v — 0, it gives the distribution of the analytic signal in 
the time domain 

M(t,l,0) = |/(t)| 2 (10) 

and for the case \i — 0, v — 1, the distribution of the analytic signal in the 
frequency domain 

M(cu,0,l) = |/M| 2 (11) 

The family of marginal distributions M(s,[i,v) contains complete infor- 
mation on the analytical signal. This may be shown directly. However it 
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M(s,n, v) 



2vrM 



exp 



i\xt 2 



its 



f{t)dt 



is more interesting to point out that there is an invertible transformation 
connecting M(s, /i, v) to the Wigner-Ville quasidistribution, namely 

M (s, fi,u) = J exp [-ik(s -lit- uco)} W(t, u) (12) 



and 



W(t, uS) = — [ M (s, fi, u) exp [—% ([it + vuj — s)\ dfi duds (13) 

Therefore, because the WV quasidistribution has complete information, in 
the sense of Eq.([|), so has M(s, /i, v). 



3 Examples 

We compute the NCT transform M(s, fi, v) for some analytic signals: 
(i) A complex Gaussian signal 



f(t) 

It has the properties 



ir) 



i/i 



exp 



a 



— r + 1 - r + % u Q t 



r 



(t) = 0, (u) = u 

(t 2 ) ~ (t) 2 = 
<^ 2 > - M a = 

2- 1 (tuj+wt)~(t) (w) 



J_ 

a 2 +/3 2 
2a 
j 



(14) 

(15) 
(16) 



This signal minimizes the Robertson-Schrodinger uncertainty relation 

•jj t > 



4 1 - r 5 



(17) 



In quantum mechanics, it corresponds to a correlated coherent state |TD| [III 
The NCT transform is 



M (s, /I, v) 



2trx? 



exp 



s - s 



\2l 



2 at 



(18) 
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with parameters 



For the case of // 







s = 


UJ U 


cos 8 

rp j Is 


sin j 



i- |i/(a-i/3) 



(19) 



evolves along the axis, changing its maximum and width 



s 



J_ 

2a 
U 



sm I 
Q 
sin 6* 



(a — i/3) — i ( - 



(20) 



Thus, we have squeezing in the quadrature components and their correlation. 
In the case /3 = 0, one has a purely squeezed state fl2| (T^j, which minimizes 
the Heisenberg uncertainty relation 



^22 
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(ii) A normalized superposition of two Gaussian signals 
f(t) = N s {A 1 f 1 (t) + A 2 f 2 (t)} 



where is 



and 



f l (t) = N i exp 



Ni 



-a^ + bit 



1,2 



"ai + a*" 


1/4 




exp 


7T 





1 (6i + 



* \ 2 



aj + a* 



(21) 
(22) 

(23) 
(24) 



The superposition coefficients being complex numbers, the normalization 
constant N* reads 



N„ 



Ai\ 2 + |A 2 | 2 + 2Re 



A,A\ h{t)f 2 {t)dt 



-1/2 



(25) 



Computing the marginal distribution M (s, /i, v) by Eq.(^j) we arrive at a 
combination of three Gaussian terms 



M(s,n,u) = N*{\A 1 \ 2 M 1 (s,fi,v) + \A 2 \ 2 M 2 (s,ii,v) 
+2Re[A 1 A* 2 M 12 {s, /i, u)}} 



(26) 
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where we have the contribution of two real Gaussian terms 

1 



Mi (s, n, v) 



cxp 



(s - Si) 2 



2a? 



i = l,2, 



and the superposition of two complex Gaussians 



M12 (s, /i, v) 



"12 



cxp 



(s - si 2 ) 2 



2a 2 
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The parameters of the real Gaussians are the dispersion 

2 



^ 2 = 2 



I/Oi - f 



Oj + Oj 



and mean 



Im Re (^) 



Re a, 



The parameters of the complex Gaussian are 



and 



^i 2 2 = 2^ 2 



s 12 



Oi + 2 



ai + 
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and the complex amplitude nyi of the complex Gaussian contribution 



Cl2 

n i2 = 7= ; exp 



27r|i/| 

(hi) Finite-time signals 
Here we consider signals 



6? 



+ 



L*2 
w 2 



4 I „, _ Mi n * _|_ Mi 



+ 



'12 



2a 2 



12 



/i(t) = iv ie - a< * 2+M , t 2i <t<t u 



which vanish for all other times and compute the NCT for one signal and for 
the superposition of two such signals. The parameters a; and 6j are complex 
numbers. The normalization constant is 



Mi = Ja,i + a* exp 
— erfc 



[ (b t + b*) 2 - 






4 [fli + a*) 




2 



'at + a* 
where erfc is the function 



tu — 



2 {en + a*) 



erfc ( yai + a* 

1/2 



2 1 



h + b* 
2 (a* + af 



erfc (y) 



e x dx 



(35) 



(36) 



Using Eq.(|5D, we arrive at the following marginal distribution 

A/? 



Mi (s, //, v) 



erfc 



-erfc 




^2i — 



vbi — is 



IjJL 



tli — 



IS 



2va>i 



(37) 



In the limit tu — > — oo, t 2i — > oo, the marginal distributions (|3T|) reduce 
to the Gaussian distribution given by fl27l) . In the case a« = 0, bi = iu>i, 
the distribution fl3"7|) describes a sinusoidal signal of finite duration. The 
normalization constant takes the limit value 



(38) 



Ni -—> \ I o, - LU) 

For a superposition of two finite-time signals 

f(t) = N s {A 1 f 1 (t) + A 2 f 2 (t)} 

with the signals fi(t) and f 2 (t) as in (|34|), the normalization constant is given 
by Eq. (|25|) with overlap integral 

+ b*f 



tb h{t)t 2 {t)dt = M^ - ^== exp 



'2^/ai + a* 2 



| erfc ^ 
— erfc 



ai + a£ 



ai + a% 



4 (aj + a*] 
h + b* 2 
2( ai + a* 2 ) 

h + b* 2 

2 (a! + aS) 



(39) 



8 



The marginal distribution for the superposition signal has the same form 
as Eq. ( p6|) but with a changed normalization constant, the distributions 
Mi (s, n, v) and M 2 (s,/i, v) given by Eq. (|37|), and an interference term 
M 12 (s, /i, i/) 



M i2 (s,fj,,u) 




(40) 



The case a 2 = corresponds to the combination of a finite time chirp and a 
finite time sinusoidal signal shown in one of the figures below, 
(iv) Graphical illustrations 

We have plotted M (s, fi, v) for some signals. In all cases we use // and v 
as in Eq.(^) with T = 1 and Q = 10. All signals are finite time signals and 
in each case we display a three-dimensional and a contour plot. 

# Figs la,b. The signal is 



/(*) 



e -i20t _|_ gilOt 







te (0,1) 1 

t i (o, i) 



(41) 



Although the number of periods, during which is signal is on, is relatively 
small, the two contributing frequencies are clearly seen in the separating 
ridges. 

# Figs 2a,b. The signal is 

- mt t g (o, i; 



/(*) 



e 




1/ 

1 3^ 
■ 4' 4' 



(42) 
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Once again the contributions separate as 9 grows, but notice the intermediate 
interference region which is a signature of the time-sequence of the frequencies 
occurrence and of their relative phase. 
# Figs 3a,b. The signal is 

' g -i(20t+10t 2 ) + e iWt 



Contrasts the signature shapes of a chirp contribution and a regular sinu- 
soidal pulse. 

Notice that all M(s, 8) values have a probability interpretation. Therefore 
all peaks or oscillations have a direct physical meaning and, as opposed to the 
time-frequency quasidistributions, we need not worry about spurious effects. 
This is particularly important for the detection of signals in noise, as we will 
see in the next example. 

(v) Detection of noisy signals by NCT 

In Fig. 4a and 4b we have plotted a time signal f(t) and its spectral density 
\f(u)\ 2 . It is really very hard to decide, from these plots, where this signal 
might have originated from. Now we plot the NCT transform (Fig. 4c) and 
its contour plot (Fig.4d) with the normalization T = 1 and Q = 1000. It still 
looks quite complex but, among all the peaks, one may clearly see a sequence 
of small peaks connecting a time around 0.5 to a frequency around 200. 

In fact the signal was generated as a superposition of a normally dis- 
tributed random amplitude and random phase noise with a sinusoidal signal 
of the same average amplitude but operating only during the time interval 
(0.45,0.55). This means that, during the observation time, the signal to 
noise power ratio is 1/10. The signature that the signal leaves on the NCT 
transform is a manifestation of the fact that, despite its low SNR, there is 
a number of particular directions in the (t,u) plane along which detection 
happens to be more favorable. The reader may convince himself of the sound- 
ness of this interpretation by repeating the experiment with different noise 
samples and noticing that each time the coherent peaks appear at different 
locations, but the overall geometry of the ridge is the same. 

Of course, to rely on a ridge of small peaks for detection purposes only 
makes sense because the rigorous probability interpretation of M(s, 9) ren- 
ders the method immune to spurious effects. 



t e (0,1) 
t£(0,l) 



(43) 
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4 NCT for other non- commuting pairs. The 
time-scale and frequency-scale cases 

The method may also be applied to other pairs of non-commuting variables 
for which, as in the time- frequency case, there cannot be a joint probability 
density. Consider the pair time-scale, where scale is the operator 



D = ]- {too + oot) = oot + % - 

In the plane (t, D) we consider the linear combination 

cos 9 

si = fit + vD = — t + vD 



T 



The relevant characteristic function is 



CW(jfe) = (jk(j*+vD)} = j f*(t)e ik ^ t+ ^f{t)dt 



= jr(e-*t)e a $**Mf(e 1 ¥t)dt 
and the NCT transform is, as before, the Fourier transform of Cj^J(k) 



-iksi 



dk 



leading to 



i 



1 

2tt\v\ 



+ 

2 



For the pair frequency- scale, (uo,D), we obtain similarly 

cos 9 



s 2 = fJM + vD 



-co + vD 



M( 2 ) 



2it\v\ 
1 

2ir\v\ 



Lo|/H^^ ?losa;) 



f{uS) being the Fourier transform of the signal f(t). 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 
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